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Motivated by recent experiments, we analyze transport of Cooper pairs through a double- island Josephson 
qubit. At low bias in a certain range of gate voltages coherent superpositions of charge states play a crucial 
role. Analysis of the evolution of the density matrix allows us to cover a wide range of parameters, inch 
situations with degenerate levels, when dissipation strongly affects the coherent eigenstates. At high noise 
levels the so-called Zeno effect can be observed, which slows down the transport. Our analysis explains certain 
features of the I—V curves, in particular the visibility and shape of resonant peaks and lines. 

PACS: 85.25.Cp, 73.23.Hk, 74.50.+r 



Among various proposals for realization of qubits, 
solid-state devices appear particularly promising since 
they can be easily scaled up to large qubit registers 
and integrated in electronic circuits Recent experi- 
ments have demonstrated quantum coherent oscillations 
in Josephson-junction devices. However, in such devices 
due to the host of microscopic modes decoherence pro- 
cesses are more difficult to control, and understanding of 
the decoherence mechanisms requires further analysis. 
Further, improvements of the quantum measurement 
procedure are needed to allow monitoring the qubit 's 
state with little influence on the qubit's dynamics before 
the read-out. Here we analyze recent experiments 013]) 
in which the dissipative dynamics of a Josephson charge 
qubit were probed by Cooper-pair transport. This ex- 
periment provides data for understanding of the dissi- 
pation in typical superconducting charge devices, and 
its analysis is similar to that for the quantum charge 
detectors. 

We focus on the analysis of Josephson circuits in 
the charge limit, in which the typical electrostatic en- 
ergy needed to charge a superconducting island (~ 
(2e) 2 /2Cs, where Cs is the total capacitance) is higher 
than the Josephson energy which controls the charge 
tunneling. If the system is biased close to a point, where 
two charge states with lowest energies are degenerate, 
at low temperatures and operation frequencies one can 
neglect the higher charge states, and the system reduces 
to two levels (qubit). The matrix element between these 
levels is controlled by the Josephson tunneling. In the 
simplest design, a Cooper-pair box 0], the quantum 
state of this qubit can be manipulated by voltage and 
current pulses [S]. The measurement of the quantum 
state can be performed, for example, by coupling the 
qubit to a single-electron transistor (SET) and moni- 
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FIG0 The double-island system 



toring its current 0. Here we study a circuit, which 
can be described as a charge qubit inside a SET. Trans- 
port in this device probes typical time scales of the qubit 
dynamics, and its analysis may show new possibilities to 
perform the read-out. Our results explain experimen- 
tally observed features of transport (the visibility and 
shape of resonant lines and peaks) and predict new spe- 
cific behavior in a low-bias regime, in which coherent 
properties of the double-island qubit are probed. 
The circuit and its description. Following the ex- 
perimental work [2 0] we study the system shown in 
Fig. n It consists of a Josephson junction, with a rel- 
atively strong coupling Ej, connected to further super- 
conducting leads via weaker junctions, with Ej C £j. 
The transport is controlled by a bias V between the ex- 
ternal leads and gate voltages V g i , V g 2 , with gate capaci- 
tances much lower than those of the junctions, C g <C Cj. 
The transport of Cooper pairs through a similar system 
with a single island between the leads (a superconduct- 
ing SET) was studied, for instance, in Refs. El E3 • 
Transport at a finite bias implies dissipation which can 
be provided by various mechanisms. Here we focus on 
low voltages and temperatures, at which the contribu- 
tion of quasiparticles is neglig ible 0. We study the 
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influence of the electromagnetic environment, i.e., ef- 
fective impedances in the circuit. Since C' g <C Cj the 
impedance of the transport voltage circuit is expected 
to dominate dissipation (see Fig.QJ. 

In Refs. the analysis was limited to the evo- 

lution of occupations of the eigenstates (the diagonal 
entries of the density matrix in the eigenbasis of the 
non-dissipativc Hamiltonian) . The dynamics were de- 
scribed in terms of incoherent transitions between these 
states. This approach is sufficient as long as fluctua- 
tions provide only a weak perturbation (the incoherent 
rates are lower than the coherent level splittings). How- 
ever, in a system with almost degenerate eigenstates 
this approach may fail, since the system crosses over to 
the so-called Zeno regime jlOj . To illustrate this con- 
cept we consider a situation relevant for the analysis 
below: often the charge transport may be described as 
a chain of transitions between various charge configura- 
tions. Under certain conditions one link in this chain is 
a pair of degenerate charge states, with the coupling 5 
between them, coupled by incoherent transitions, with 
rate ~ T, to further states. As long as 5 ^ T transport 
within the pair is fast, and the current magnitude is set 
by r. However, if the coupling S becomes weaker than 
F, the dynamics change dramatically: frequent 'obser- 
vation' (fast dephasing) by the transitions destroys the 
coherence and slows down the evolution; the system is 
blocked for a long time in one of the charge states in 
the pair, with the typical transition rate ~ S 2 /T, which 
now sets the current magnitude. The density matrix of 
the two-state system quickly becomes diagonal in the 
charge basis, while in the eigenbasis diagonal and off- 
diagonal entries of the density matrix are strongly cou- 
pled (cf. Ref. PP). In order to describe the behavior of 
the system in both limits, we analyze the system using 
the master equation for the evolution of all entries of 
the density matrix. 

We describe the state of the system by the charges 
eni , eri2 of the central islands and introduce the charge 
em transferred across the system of three junctions (see 
below for a precise definition). In the Hamiltonian, 



H — He + Hj + i?diss , 

the charging part is given by 

_ (en_+C g U g _) 2 {en + + C g V g+ f 
C 4(3Cj + C g ) 4(Cj + C g ) 

where n± = n\ ± na, V g ± = V g i ± V g 2 and 

en_(2Cj + Cg) en + C g 



lint 



(1) 



V, (2) 



is the charge operator that couples to the voltage source. 
The Josephson part of the Hamiltonian is 



H, 



cos0 2 )-£ J cos(0 2 -0i+*m), (4) 



where 9\, 9 2 are the phase drops across the left and 
the right junctions, respectively, and exp(i^ m ) : \m) i— > 
\m + 2) is the counting ladder operator. One can see 
from Eq. J2Jl that tunneling of a Cooper pair across the 
central junction changes m by 2. 

Finally, the dissipative part of the Hamiltonian reads 
(cf. HUE]): 



■"diss 



(Qint - qf 
2Ci n t 



2C a 



2L r 



(5) 



2(3Cj + C g ) 2(Cj + C g ) 



+ em (3) 



Here <J) is the phase drop across the impedance Z(u>) 
and Ci„t is the capacitance between its leads. The lin- 
ear environment is presented here as a parallel connec- 
tion of LC-oscillators, with the constraint Z^ito) = 
J2 a ,± [ iL <x(u ± + «0)] _1 , where uj a = l/y/L a C a . 

We obtained a Hamiltonian description in terms of 
the phases 6i, 62, 4> an d the conjugate charges eni, 
— eri2 and q, the latter being the total charge passed 
through the voltage source relative to the equilibrium 
charge C\ nt V on the plates of the capacitor Ci n t- At the 
relevant low frequencies the interaction with the bath 
reduces to Hi nt — — Q m t5V, where 8V = (q — m)/Ci nt 
is the fluctuating part of the transport voltage. The 
dissipative Hamiltonian of the form (jjjj provides for 
the proper high-frequency regularization of the effective 
bosonic bath, with cut-off frequency u> c — (RC- m t) , 
where R is the real part of Z(oj). 

Qualitative analysis of the low-voltage reso- 
nances. In this section we provide qualitative analysis 
of the transport properties and illustrate the discussion 
by the results of numerical simulation described below. 
We study resonances at transport voltages below the su- 
perconducting gap, eV < A but assume that the voltage 
is high enough so that the features related to the super- 
current through the system are not relevant. The dis- 
cussion and figures correspond to a positive bias V > 0. 
To understand the origin of possible resonances, let us 
first discuss the stability diagram for the charge states 
(see Fig. EJ, neglecting the Josephson couplings. 

In the unbiased case the stability conditions define 
a honeycomb pattern in the gate-voltage plane. Inside 
each hexagon a certain charge state has the minimal 
energy. At the vertices three charge states are degen- 
erate. When a transport voltage V is applied these 
points grow into triangles, within which the system is 
unstable w.r.t. sequential tunneling of Cooper pairs: 
|0,2,m) -> |0,0,to) -> \2,0,m) -> |0,2,m + 2)... [In 
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and coherently through two other junctions. The co- 
herent couplings Ej exceed the rate of incoherent tran- 
sitions, which can be evaluated using the golden rule: 



FIG|5| The honeycomb stability diagram of charge 
states. The solid dots and dashed lines denote the res- 
onance peaks and lines, respectively 
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FIG[3] Three charge levels at resonance. The (00) and 
dotted arrows denote a passage along a resonant line 



the experimentally relevant limit of Josephson couplings 
and temperatures below the charging energy, in the 
vicinity of one vertex only three charge states |ni,n 2 ) 
are relevant: |2, 0), |0, 2) and |0, 0).] However, this gives 
a low current since the incoherent tunneling through the 
left and right junctions Ej is slow. 

A much higher current can be achieved in resonant 
situations. One can expect resonant points (peaks) and 
lines in the Vg±-plane. At the peaks, defined by two 
constraints on V g ±, three charge states are in resonance. 
On the lines, only two charge states are degenerate. The 
resonant conditions determine the positions of possible 
peaks and lines. To evaluate the current at the res- 
onant peaks, we note that for typical parameters the 
bottleneck of transport is associated with the incoher- 
ent transitions between triples of resonant states; the 
rate of these transitions is given by the golden rule and 
defines the current. However, the analysis of the shape 
of the peaks/lines (the decay of current away from res- 
onances) is more subtle. It may require the analysis 
of the Zeno regime and of the crossover to this regime. 
Below we develop a suitable master-equation approach. 
We begin with a qualitative discussion of the results. 

Consider, for instance, the three-state resonance 
shown in Fig. |3 which corresponds to the upper ver- 
tex of the triangle in Fig[21 In this case the Cooper 
pairs tunnel incoherently in the central junction only, 



r r 



4tt R El 
Y R Q 2eV 



(6) 



where Rq = h/(2e) 2 and we assumed T <C 2eV. This 
rate defines the current magnitude at resonance, 2er r . 

Tuning the gates away from this resonance peak, 
one may still keep two levels degenerate along a reso- 
nant line. For instance, one may lift the state |0,0,m) 
with respect to the degenerate |0, 2, m) and |2, 0, to) (see 
FigljjJ if |0,0) descends, the system may get Coulomb- 
blocked in this state). In this configuration the trans- 
port involves a second-order coherent tunneling (cotun- 
neling) |0, 2, to) — > |2,0,to) and incoherent relaxation 
|2,0, to) — > |0, 2, to + 2). To estimate the current, we 
evaluate the second-order coherent coupling between 
|0,2,to) and \2,0,m) and find S ~ E^/AE, where AE 
denotes the distance to the |0,0)-state (see Fig|3J). As 
discussed above on p|2 the current remains to be 2er r 
as long as 6 > T r . However, for S < T r the system is in 
the Zeno regime, and the relaxation rate |2,0) — > |0, 2) 
defines the current ~ 2e8 2 /T T . Thus along the resonant 
line the current stays at the peak level and then drops 
fast. The deviation from the peak at which the current 
drops can be estimated from the condition S ~ T r ; fur- 
ther behavior is governed by the Zeno physics. For the 
typical parameters (see below) this gives a very short 
line (it is also very narrow, cf. below). This may explain 
why this resonant line was not detected. 

If the threefold degeneracy in Figgis lifted in other 
ways (with two states still in resonance) the transport 
involves higher-order incoherent processes and the cur- 
rent is much weaker However, there exist other res- 
onant peaks, which are located at two lower vertices 
of the triangle: one can say that in FigEl the volt- 
age drops at the central junction, but it can also drop 
at the left/right junction. The respective rate of in- 
coherent Cooper-pair tunneling can be evaluated using 
Eq. © with the substitution Ej — > Ej, i.e., the current 
at these peaks is much lower. However, our analysis 
shows that the resonance lines originating from these 
peaks are much longer (and wider, cf. below) and may 
reach the neighboring hexagons' vertices, as it was in- 
deed found experimentally. The reason is that at these 
peaks the incoherent rate is much lower, the coherent 
coupling stronger, and it takes a longer distance away 
from the peak for the coherent coupling to fall below the 
incoherent rate (crossover to the Zeno regime). Thus we 
find, in agreement with experiment, that only oblique 
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FIG|1] Resonances with the double- island's eigenstates. 
The right panel shows the peaks' positions in the V s ±- 
plane. A peak emerges at V = Ej/2e (the solid dot 
in the middle) and splits as the bias V increases. The 
dashed lines show the cotunneling resonances 

(but not horizontal) resonant lines should be visible and 
allows us to evaluate the shape of the resonances. 

The widths of the resonant lines were evaluated in a 
similar way, with results in, at least, semi-quantitative 
agreement with experiment. We remark that the width 
is not set by the condition of resonance as such (which 
requires a charge-level splitting lower than the coupling 
and would define very narrow lines |3]). In fact, dur- 
ing the separation of two resonant states the transport 
changes from coherent to incoherent. At this crossover 
point the incoherent rate is higher than the respective 
T r . Only at a longer distance from the line it drops be- 
low F r and slows the transport. The respective width 
scales linearly with V, similar to the experiment [2]- 

So far we analyzed transport at voltages V much 
higher than the Josephson couplings and used the 
charge basis. Now we focus on transport at lower volt- 
ages and find that due to the coherent Josephson cou- 
pling of the charge states the triple resonance of FigEl 
appears only at voltages above a certain threshold. 

At lower 2eV ~ Ej it is convenient to work in the 
eigenbasis of the double island. Near the triangle in 
Fig.[21thc difference U in charging energies of the states 
|2,0,m) and |0, 2, m+2) is small, and one finds the 
ground and excited eigenstates of the double island, 

\g,m + 2) =cos7|2,0,m) + sin 7 |0, 2, m + 2) , 
|e, to) = — sin7 |2, 0, m) + COS7 |0, 2, to + 2) , 

where tan 27 = Ej/U. 

The respective resonance configuration is shown in 
FigEl Since the minimal energy splitting between the 
ground and excited states is Ej, the resonant condi- 
tions of Fig. H require 2eV > Ej. At V = Ej/2e 
the peak is located at the lower side of the trian- 
gle (see Fig. 0J. Above this threshold the equation 
E c - E g = £j/sin27 = 2eV (0 < 7 < tt/2) has two 
solutions, and the peak splits: The main peak with 
7 > 7r/4 enters the triangle, and the other, secondary 
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FIGEl I(Vg+,V g -) for V = 23/iV and T = 50mK 



peak (7 < 7r/4) leaves it. At strong bias 2eV ^> Ej 
the main peak reaches the upper vertex of the triangle, 
while the secondary becomes very narrow and joins one 
of the oblique resonant lines. (Notice that the triangle 
itself slides and grows with the increase of V .) 

Let us estimate the current magnitude at these 
resonances. The relaxation rate |e, to) — > |<7, to + 2) 
is given by Eq.©, and the matrix element between 
the states |e,m) and \g,m) to |0, 0,m) due to Hj is 
E'coupi = (Ej/2) sin7, of order Ej for the main reso- 
nance • and weaker, ~ EjEj /(2eV), for the other one. 

If Ej (R/Rq) Ej the incoherent relaxation inside 
the double-island is the bottleneck (the slowest stage) 
of the transport for both peaks, T r -C -Ecoupb i-e., the 
peak height is / ma x ~ 2er r . However, the peaks' sizes 
are different due to different -Ecoupi and can be found 
from the analysis similar as above. The external peak o 
is much narrower at 2eV 3> Ej . 

For Ej < (R/Rq) Ej one finds that T r > E coup i for 
the secondary peak, and also for the main resonance at 
voltages V close to Ej/2e. The Zeno effect is expected 
under these circumstances ^0l> the transport is slowed 
down, and transitions in the outer junctions define the 
current J max ~ 2eJ5^ oupl /r r . 

Master equation and numerics. The dynamics re- 
duce to propagation along the chain of eigenstates with 
decreasing energy and growing to. To evaluate the cur- 
rent we analyze the dynamics of the reduced density ma- 
trix <t, retaining the indices 711,712, to and tracing over 
the environment's degrees of freedom. Using the real- 
time Keldysh diagrammatic technique (cf. Ref. |lHU14j ). 
we find the master equation 



~-&{t) - L &{t) = / dt' m-t')a(t'). 



(7) 
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with the bare Liouville operator Lq = i[-,Ho], Hq = 
Hq + Hj. In the first (Born) approximation we obtain 

E(t) = a'{t)L int e Lot L int - ia" (t)L int e Lot M int , (8) 

where a(t) — a'(t) + ia"(t) is given by 

a(t) = (2ef(6V(t)5V(0)) = / ^ jjgg^ , 0) 

the low-frequency spectral density J(w) = 2t:lu R/ Rq, 
and L int = i [•> Qint/2e], M int = i [•, <9 in t/2e] + . The 
last term in Eq. © violates the translational symme- 
try m — » m + 2. The invariance is restored after a 
regularization, due to the counterterm Q? nt /2C ln t in 
Eq. © (cf. Ref. E3). 

We label the entries of the self-energy matrix X by 
four triples and z/^, where e.g. 2/~ = in{ , , to - ). 
Here the sign qp refers to a Keldysh branch; the 
(un)primed indices refer to the time t' (resp. t). Most 
of these indices vary over finite ranges. Indeed, only 
the lowest charge states m, n-i participate in the low- 
frequency dynamics, and strongly off-diagonal entries, 
with large m~— m + and m'~— m' + , are suppressed. The 
regularized self-energy is translationally invariant and 
does not depend on the sum m~+ m + + m' - + to' + . The 
Fourier transform w.r.t. (m~+ to + — to' - — m l+ )/2 gives 
a finite matrix for each value of k. 

We use the Laplace-transformed master equation: 
sa(k,s) — (To(fc) = II(s, k)a(k, s) to find the current 
I = s 2 (m(s))\ s ^o, where (m(s)) = ie\.Tr a(k = 0, s) = 
iTr (s - iTr 1 ^!! (s - n)- 1 & U-o- Here II(fc, s) ee 
Lo(k) + s) and <To is the initial condition. The nu- 
merical analysis can be simplified by taking the needed 
derivatives analytically and working in the eigenbasis 
of Hq. We ascribe a counting index m to eigenstates 
(rather than charge states) and organize them into zones 
with fixed values of to |S] . The eigenstates of the to- 
tal Hamiltonian Q have only a finite TO-spread, and one 
can use rh to evaluate the dc current. 

In our analysis we used the following parameters: 
Cj = 0.8 fF, C g = 8aF, Ej = 25 mK, Ej = 0.5 K, 
R ~ 50il 3 . Fig[5] shows the shape of two peaks and 
resonance lines for the bias V just above the threshold 
Ej/2e ~ 21.5/iV, in agreement with estimates above. 
Discussion. In our calculation we neglected the influ- 
ence of the 1// noise due to background-charge fluc- 
tuations. This very-low-frequency noise dominates the 
pure dephasing (that leads to energy fluctuations with- 
out transitions; its Ohmic part is included in our nu- 
merical analysis). Our estimates show that these effects 
should not change the results substantially. 

In conclusion, using the methods that allow to 
cover the (Zeno) dynamics of coherent systems under 



strong dissipation, we analyzed the Cooper-pair trans- 
port through a double-island structure. We find sep- 
arate peaks and resonant lines, whose visibility and 
shapes match the experimental observations. We fur- 
ther predict a double-peak structure near a thresh- 
old transport voltage, observation of which would be a 
probe of coherent properties of the double- island qubit. 
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